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Abstract. This contribution compares several different approaches allowing one to derive macroscopic 
trafhc equation directly from microscopic car-following models. While it is shown that some conventional 
approaches lead to theoretical problems, it is proposed to use a smooth particle hydrodynamic approach and 
to avoid gradient expansions. The derivation circumvents approximations and, therefore, demonstrates the 
large range of validity of macroscopic traffic equations, without the need of averaging over many vehicles. It 
also gives an expression for the "traffic pressure", which generalizes previously used formulas. Furthermore, 
the method avoids theoretical inconsistencies of macroscopic traffic models, which have been criticized in 
the past by Daganzo and others. 
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1 Introduction 

In order to describe the dynamics of traffic flows, a large 
number of mathematical models has been developed. The 
analysis of the spatio-temporal features and statistics of 
traffic patterns has often been done with methods from 
statistical physics and non-linear dynamics. An overview 
of modeling approaches and methods is, for example, given 
in Refs. [UElElll] , among them cellular automata, "mi- 
croscopic" car-following models, "mesoscopic" gas-kinetic, 
and macroscopic traffic models. 

Cellular automata can often be interpreted as discretized 
versions of car-following models, while gas-kinetic models 
have frequently been used to derive macroscopic from mi- 
croscopic models. Such derivations were driven by the de- 
sire to improve phenomenological specifications of macro- 
scopic traffic models [5l[6j[7] , which were criticized to have 
unrealistic properties However, the derivation of gas- 
kinetic models from car-following models usually simplifies 
the interactions among vehicles by a coUisional approach 
assuming immediate braking maneuvers. Moreover, the 
derivation of macroscopic traffic models from gas-kinetic 
ones terminates an infinite and poorly converging series 
expansion, which replaces dynamical equations for higher 
moments of the velocity distribution by simplified equilib- 
rium relationships [8j. 

Although this leads to macroscopic equations which 
work well in most theoretical and practical aspects ,9 , 
the implications of the approximations are hardly known. 
Moreover, the approach seems to require an averaging over 
at least 100 vehicles for each speed class and spatial loca- 



tion. While this constitutes no problem for gases with 10^^ 
particles within a small volume, for traffic flows this would 
require an averaging over spatial intervals much greater 
than the scale on which traffic flow changes. Hence, it is 
not well understood, whether or why macroscopic traffic 
equations can be used at all. 

In this paper, we will therefore focus on attempts to 
derive macroscopic traffic equations directly from micro- 
scopic ones. Doing so, we will compare three different ap- 
proaches: First, we study the gradient expansion approach 
in Sec. [21 Second, we turn to the linear interpolation ap- 
proach in Sec. [21 Third, we discuss the smooth particle hy- 
drodynamics approach in Sec. [H and compare the results 
with macroscopic traffic models such as the Payne model, 
the Aw-Rascle model, and a non-local traffic model. In 
the Conclusions, we summarize and discuss our results, in 
particular with regard to the mathematical form of the 
traffic pressure and the theoretical consistency of macro- 
scopic traffic models. 



2 The Gradient Expansion Approach 



Already in the 1970's, Payne [TUlfTTj used a gradient ex- 
pansion approach to derive a macroscopic velocity equa- 
tion complementing the continuity equation 



dt dx 



p{x,t)V{x,t) 



(1) 
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It relates the vehicle density p(x, t) at location x and time 
t with the average velocity V{x,t) or the vehicle flow 

Q{x,t)=p{x,t)V{x,t), (2) 

respectively, and describes the conservation of the number 
of vehicles [12] ■ 

Payne derived his model from Newell's car-following 
model [13] 

v,{t + t) ^ v,{d^{t)) , (3) 

which assumes that the speed Vi{t) of vehicle i at time t 
will be adjusted with a delay of r to some optimal speed 
Vo, which depends on the distance di{t) — Xi-i(t) — Xi{t) 
between the location of the leading vehicle Xi-i (t) and the 
location Xi{t) of the following car. 

Payne identified microscopic and macroscopic veloci- 
ties as follows: 



term, as it reflects the reaction of drivers to the traffic situ- 
ation in front of them. The relaxation term [Vc{p) — V]/ At 
delineates the adaptation of the average velocity V{x,t) 
to the density-dependent equilibrium velocity Vc (p) with a 
delay t. 

Other authors have applied similar gradient expan- 
sions to the optimal velocity model defined by 



dvi{t) 
dt 



Vo{di{t)) - Vi{t) 



(10) 



Vr{t- 



V{x 

V{x,t) + VT 



VT,t + T) 

dV{x,t) 



dx 



dV(x,t) 
dt 



(4) 



Then, Taylor approximations (gradient expansions) were 
used in several places. For example, Payne substituted the 
inverse of the distance di to the leading vehicle by the 
density p at the place x-t-dj [t) /2 in the middle between the 
leading and the following vehicle. In this way, he obtained 



with ddi/dt = Vi-i{t) — Vi{t), see e.g. Refs. fl5','16'. Equa- 
tion pil| results from the Newell model ([3]) by a first-order 
Taylor approximation Vi (t + r) « Vi (t) + r dvi / dt. Regard- 
ing the derivation of macroscopic traffic equations from 
the optimal velocity model, it is also worth reading Refs. 

One weakness of the gradient expansion approach is 
that it implicitly assumes small gradients in order to be 
mathematically valid. It is well-known, however, that many 
microscopic and macroscopic traffic equations give rise to 
emergent traffic jams, which are related with steep gradi- 
ents. That would require the consideration of higher-order 
terms and lead to macroscopic traffic equations that are 
not anymore simple and well tractable (even numerically). 
Let us, therefore, study other approaches to determine 
macroscopic from microscopic equations. 



d^it) 



d,it) 



1 

2^' 



3 The Linear Interpolation Approach 



p{x,t) + 



1 dp{x,t) 



(5) 



2p dx 

When defining the so-called equilibrium velocity Vc{p) through 



V.ip) - 



or 



V, 



^ ] = Vo{d^) , 



(6) 



The optimal velocity model may be also written in the 
form 

-^=a,{t)^ -l^ + f{d,{t)), (11) 

where ai{t) denotes the acceleration, the "desired ve- 
locity" or "free speed" , and 



a first order Taylor approximation and Eq. ([5]) imply 
1 



Vo{d,{t)) = V, 



m) = ^^^^^^^ < 



(12) 



d,{t) 
Vo{p{x,t)) 



1 



dV,{p) dp{x,t) 



2p(x, t) dp 



dx 



(7) 



Starting from the previous equations, one finally arrives 
at Payne's macroscopic velocity equation 



the repulsive interaction among the leading vehicle z — 1 
and its follower i. 

In Ref . [17] , it has been suggested to establish a micro- 
macro link between microscopic and macroscopic traffic 
variables by the definitions 



1 



ay ^dv_ _ 1 

dt dx T 



D{p) dp 



V{x,t) 



p dx 

where we have introduced the density-dependent diffusion 

dV,{p) 



(8) p{x,t) 



Xi{t) - Xi+l{t) 



[Xi-l{t) 



Xi 

1 



Xi^l{t) - Xi{t) 



i{t)-x.,{t) 

' X^{t)\ 



D{p) 



1 dV,{p) 

2 dp 



dp 



> 0. 



(9) 



The single terms of Eq. ([8]) have the following interpreta- 
tion: The term VdV/dx is called the transport term and 
describes a motion of the velocity profile with the vehi- 
cles. The term — [D{p)/{pAt)]dp/dx is called anticipation 



V{x,t) 
A{x,t) 



Xi-l{t) - Xi{t) 

■H{t)[xi-i{t) -x] + v.,^i{t)[x - Xi{t)\ 



ai{t) [a 



Xi^l{t) 

-i{t)-x\ 4 



Xt{t) 
1i-l(t)[x - 



Xi^l{t) - Xi{t) 



(13) 
(14) 
(15) 
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These definitions assume that the macroscopic variables in 
the vehicle locations x — Xi(t) would be given by the mi- 
croscopic ones, while in locations x between two vehicles, 
they would be defined by linear interpolation. 

Let us consider the consequences of such an approach. 
For this, we determine the partial derivative of 

r<f .^9zit)[xt-iit)-x]+g,-i{t)[x-x^{t)] 

x.^.it) - x.it) (^^^ 
with respect to x, which gives 

dGix,t) ^ -g^(t) + g,_,{t) ^^^^ 

dx Xi-i{t) - Xi{t) 

for any specification of gi{t)^ for example, gi{t) = Vi(t). 
The partial derivative with respect to time is 

dGix,t) _ ^[x,_,it)-x] +g.it)'-^ 



dt 



dt 



Xi-^l{t) - Xi{t) 

X ~ Xi{t)] - gi-i{t) 



dxjjt) 
dt 



Xi-l{t) - Xi{t) 

('^-^)<>.Wh-iW-^l 



( 



dxi-i [t) 
dt 



[Xi-i{t) - X^(t)\^ 

gi-i(t)[x 



dxi(t) 
dt 



(18) 

[Xi-i{t) - Xi{t)\ 

For gi{t) — Vi{t) = dxi/dt and with dvi/dt = ai{t), this 
formula simplifies to the following expression: 

dV{x,t) _ a,{t)[x,^i{t) 

dt Xi^i{t) ~ Xi{t) 

at-i{t)[x - Xi{t)\ - Vi^i{t)vi(t) 

Xi-i{t) - Xi{t) 
_ Vi-l(t) - Vi(t) 
Xi^l{t) - Xi{t) 

^Vi{t)[xi^l{t) -x\ + Vi^l{t)[x - Xi{t)\ 
Xi-l(t) - Xi{t) 

^Aix,t)^^^^^Vix,t). (19) 
As a consequence, we find the exact relationship 

^ + V(.,,)^..4(,M)^ (20) 

This would be fully compatible with Payne's macroscopic 
traffic equation 



if 



1 



A{x,t) = - V,{p)-V{x,t) 



Djp) dp 
Tp{x,t) dx 



(21) 



However, the expression for (7i(t) = l/[xi-i{t) — Xi{t)] does 
not simplify in a way that would finally lead to the con- 
tinuity equation ([1]). Therefore, a micro-macro link based 
on the linear interpolation (jl6[) of the microscopic vari- 
ables gi (t) does not exactly imply the conservation of the 
number of vehicles, i.e. it is theoretically not consistent. 
Nevertheless, it works surprisingly well in practise |17j . 



4 The Smooth Particle Hydrodynamics 
Approach 

4.1 Derivation of the Continuity Equation 

In this section, we will derive macroscopic traffic equa- 
tions directly from microscopic ones. We will start with 
the derivation of the continuity equation from the equa- 
tion of motion dxi/dt = Vi, using a "trick" that I learned 
from Isaac Goldhirsch. For this, we represent the loca- 
tion Xi{t) of an element i in space by a delta function 
5{x — Xi{t)), which may be treated here like a very narrow 
Gaussian distribution. Moreover, we introduce a symmet- 
rical smoothing function 



s{x' — x) — s{\x' — a;|) = s{x — x') , 



(22) 



for example, a Gaussian distribution with a finite variance 
or a differentiable approximation of a triangular function 
or a rectangular one. The smoothing function shall be nor- 
malized by demanding 



dx' s{x' - x) = 1 (23) 
for any value of x. With this, we define the local density 

oo 

p{x,t)= J dx' six' - x)Y^S{x' - x^{t)) (24) 

-oo * 

= "^s{x^{t)-x). (25) 

i 

Herein, we sum up over all particles i. Note that the re- 
placement of the conventional formula d{xi{t) — x) for 
the vehicle density by the formula s{xi(t) — x) corre- 
sponds to a substitution of point-like particles by "fuzzy" 
particles, which is the idea behind smooth particle hydro- 
dynamics. 

Now, we define the average velocity V{x,t) as usual 
via a weighted average with the weight function S{x' — 
Xi{t))s{x' — x): 



V{x,t) = 



oc 

/ dx' ^v.i{t)6{x' - x,{t))s{x' - x) 
i 

oo 

/ dx' y~^(5(a; — Xi{t))s{x' — x) 

-oo ^ 

oo 

/ dx' ^Vi{t)5{x' - Xi{t))s{x' - x) 
p{x,t) 

^^Vi(t)s{Xi{t) - x) 

i 

^s(x,(t) -x) 

i 

i 

p{x,t) 



(26) 
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This implies the well-known fluid-dynamic flow relation- 
ship 

Q{x,t) ^ p{x,t)V{x,t). (27) 



Differentiation of Eq. ([24|) with respect to time and appli- 
cation of the chain rule gives 



dpjx, t) 
dt 

OO 



dxi 
~dt 



d_ 

dx 



-;5{x' ~ Xi{t) 



OO 

j dx' J2^m6(^ x' - Xi{t) 



s{x' — x) 
, (28) 



where we have applied partial integration to obtain the 
last results. That is, we have used the theorem 



dx' 



u{x)v(x) 



v{x') 
u{x') 



9 , , ' 



(29) 



considering the vanishing of the flrst term after the equal- 
ity sign due to the vanishing of u{x)v{x) at the bound- 
aries. Taking into account the symmetry of the smooth- 
ing function s{x' — x), we may replace ds{x' — x)/dx' by 
—ds{x' — x)/dx, which finally yields Eq. ([T]) as follows: 



dp{x,t) 
dt 



d_ 

dx 

d_ 

dx 



OO 

j dx' Vi{t)5(^x' — Xi{t))s{x' — x) 



pix,t)V{x,t) 



(30) 



To obtain this desired result, we have finally applied the 
definition (^5)1 of the average velocity V{x,t). As a con- 
sequence of this, the validity of the continuity equation 
does not require an averaging over large numbers of enti- 
ties, i.e. macroscopic volumes to average over. This makes 
the equation absolutely fundamental and explains its large 
range of validity. 



4.2 Derivation of the Macroscopic Velocity Equation 

In order to derive the equation for the average velocity, 
we start by deriving the formula 



P 



ix,t)Vix,t) ^J2v,{t)s{x,it) 



(31) 



for the vehicle flow with respect to time. This gives 
^^[pia:,t)V{x,t)]=Y^^4x.{t)-x) 



+ E^'(*)a^^(^^(^)- 

i 
i 

- ^E[^*(^)]'^(^^(^) 



dxj {t) 
dt 



(32) 



Introducing Svi{x,t) = Vi(t) — V{x,t) and defining the 
velocity variance 



e{x,t) 



J dx' J2^[y^{t)~V{x,t)]■'S{x' -X,{t))s{x' -X) 

— OO 

OO 

/ dx' Y.^S{x' ~X,{t))s{x' -X) 

— OO 

j:^[v,{t)-Vix,t)]^six,it)-x) 
^^s(a;,(t) - x) 

J2^[Sv^{x,t)]'^s{x^{t) - x) 

p{x,t) 



(33) 



similarly to the average velocity (|26p . we can make the 
decomposition 

J2[vM^s{x,it)~x) 

i 

= ^[y(x, t) + 6v,ix, t)fs{x,it) - x) 

i 

= ^{[F(x,i)]2 + 2F(a:,t),5w,(a:,t) 

i 

+ {bv.,{x,t)f'^s{x^(t) - x) 
= p{x, t){V(x, t)f ^ 2p{x, t)V{x, t) [V{x, t) - V{x, t)] 
+pix,t)e{x,t), (34) 
where we have considered 

Svi{x, t)s(xi{t) — x) 

i 

= E h(*) - ^(^' *)] K^»(^) - ^) 

i 

^ Q{x,t) ~ p{x,t)V{x,t) =0, (35) 
see Eqs. (^5)) and (^5]). Altogether, we get 

^ [p{x, t)V{x, t)] = ~^{p{x, t) [V{x, tf + 6{x, t)] } 



dx 

^a^{t)s(xi{t) - x) 



(36) 



Now, we carry out the partial differentiation applying the 
product rule of Calculus. Taking into account 

P(--0^ = -Vix^tf-^ + i[pix,t)Vix.t)] 



dt 



dt dt^ 



(37) 
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and 



d_ 

dx 



[[p{x,t)V{x,t)]V{x,t)] 
dV 



p{x,t)V{x,t) 
d 



dx 



+ y{x,t)Q^ p{x,t)v{x,t) 



(38) 



we obtain with Eq. ([36 
.dV{x,t) 



p{x,t)- 
-V{x,t) 



dt 

dp{x,t) 



d 



V{x,t)—[p{x,t)V{x,t)\ 



- p{x,t)V{x,t) 
+ ^a,{t)s{xi{t) 



dt ' dx 

dV{x,t) d 



dx 



(39) 



Inserting the continuity equation ((30|) for dp/dt and divid- 
ing the above equation by p(x, t) finally gives the velocity 
equation 

+ V[x, t)- 



dt 
1 



dx 



d 



p{x, t) dx 



[p{x,t)e{x,t)\ 



Inserting Eq. pT|) for ai{t), we find 
ai{t)s{xi{t) - x) 



(40) 



E 



s{xi{t) - x) 



v° - V{x,t) 



Y.f{dm)s{x,{t)-x). (41) 



For further simplification, let us now specify the smooth- 
ing function by the rectangular function 

^ * ' 2 1 otherwise, ^ ' 

with a large enough smoothing window of length Ax = 
2/ Q. Then, the number of vehicles i within the smoothing 
interval [x — 1/ q,x + 1/ g] is expected to be p Ax = 2p/ g, 
where p represents the average vehicle density in this in- 
terval. Therefore, 



2pg 
g 2 



(43) 



which shows the consistency of this approach. 



If the smoothing parameter g is specified via the in- 
verse vehicle distance 



1 



1 



Q = Qk = -r = 



Xk 



= p{x,t) for Xk<x<Xk-i, 

(44) 

the smoothing window of length Ax — 2/ g will usually 
contain only two vehicles k—1 and k with Xk < x < Xk-i- 
With this, the sum over i reduces to two terms with i — k 
and i — k — 1 only. This finally yields 

V{x,t) ^J2v,{t)s{x,{t)~x) 

i 

= Vk{t)s{xk{t) - x) + Vk-l{t)s{xk-l{t) - x) 



^[Vk^l{t)+Vk{t)\ 



p{x,t) 



Ufe_l(t) + Vk{t) 



(45) 



and, considering Eq. (pi)) . 



2 \9kJ 2 \gk^i 
^ p{x,t) ^ ( 1 



p{x,i) 



p{x,t) 



f 



p{x + l/p, t) 



(46) 



In summary, the macroscopic velocity equation correspond- 
ing to the optimal velocity model corresponds tcQ 



dt 

- 1 d 
p{x, t) dx 

+ 2^(k^ 



dx 



[p{x,t)9{x,t)\ 



v° - V{x,t) 



1 



^ 2^ \p{x + l/p,t) 



(47) 



It should be noted that this equation is non-local due to 
the dependence on x + 1/ p{x, t). This reflects the anticipa- 
tory behavior of drivers, who react to the traffic situation 
ahead of them. From the point of view of traffic simulation, 
the non-locality does not constitute a problem. Non-local 
traffic models such as the gas-kinetic based traffic model 
summarized in Appendix |X] can be even numerically more 
efficient than local ones with diffusion terms, that would 
result from a gradient expansion. In fact, the reason for 
the numerical inefficiency of explicit solvers for partial dif- 
ferential equations is the diffusion instability, which must 
be avoided by small time discretizations llSj. As pointed 



^ If another smoothing function is applied, the last term of 
Eq. (|47p is replaced by a similar weighted mean value, as Eq. 
(|4ip reveals, but the essence stays the same. That is, the way 
of looking at the microscopic equations (i.e. the way of defin- 
ing the density and velocity moments) potentially has some 
influence on the dynamics, but it is expected to be small. 
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out by Daganzo [19J, a diffusion term also implies theo- 
retical inconsistencies such as the occurence of negative 
velocities at the end of jam fronts. Therefore, it should be 
underlined that numerical inefhciencies and theoretical in- 
consistencies can be avoided by working with the non-local 
velocity equation rather than with the gradient expansion 
of it. which will be looked at in the next section. 



4.3 Comparison with Other Macroscopic Traffic 
IVIodels 



If the velocity variance 9 is zero, this model corresponds 
exactly to Payne's macroscopic traffic model with the pres- 
sure term p^fTT] 



P{P) = 



2t 



(54) 



It should be noted that the gradient dP/dx = [dP{p)/dp] 
■dp/dx of this pressure becomes zero, whenever the den- 
sity becomes zero or maximum, as the derivative of 
dVc{p)/dp vanishes in these situations. 



According to the discussion above, a gradient expansion 
is acceptable in case of small gradients, e.g. when a linear 
stability analysis is performed. It is also useful to compare 
different macroscopic traffic models. For this purpose, let 
us carry out a Taylor approximation of first order. It gives 



1 



^p{x + l/p,t) 
1 



/ 



p{x,t) \ dx p{x,t)^ _ 

1 \ df{d) ( dp{x,t) 1 



^p{x,t) 



dd 



dx p{x,t)^ 



(48) 



where we have applied the geometric series expansion 1/(1 — 
z) 1 + z + . . . Note that the relation p = 1/d and 

Veip) = K = Void) =V"+ Tf{d) =V"+Tf 

(49) 

imply 



dfjd) ^ / d V,ip)-v^ \ dp ^ 1 dV,{p) ^ (l_ 
dd \dp T ) dd T dp \ <P 

T dp 

Therefore, using Eq. we finally obtain: 



(50) 



^s(x.(t)-x)/(0«p(a:,t)/f-i 



^ + y(.,.)^ 
at ox 

1 dP{x,t) ^ V,{p)--V{x,t) 



p{x,t) dx 



+ 



(53) 



4.3.1 The Macroscopic Traffic IVlodel by Aw and Rascle 

Note that Daganzo has seriously criticized macroscopic 
traffic equations of the type ([S3]) [5D]. For example, he 
considered the case of a vehicle queue of maximum density 
p — pjam and speed V = Vc(pjam) = 0, the end of which 
was assumed to be at some location x = xq. Then, for 
the last vehicle in the queue, Eq. ([53]) predicts V = 
and dV/dt — dV/dt + VdV/dx < 0, i.e. the occurence 
of negative velocities, if pressure relations such as P = 
P^o ~ rjt^dV/dx with non- negative parameters and r/o 
are assumed [7], as it was common at the time when Ref. 
[T9] was published. 

In order to overcome Daganzo's criticism. Aw and Ras- 
cle [5D] have proposed the macroscopic velocity equation 

j^[V+p{p)] + V—[V+p{p)]^G (55) 

with p{p) = p'^ . Let us study, how this model relates to 
the previous macroscopic models. For this purpose, let us 
apply the chain rule of Calculus to obtain 



di 

dp{p) dp{x,t) 



dp 



dt 



dx 
V{x,t) 



dp{p) dp{x,t) 



dp 



dx 



(56) 



Inserting the continuity equation for dp/dt on the 
right-hand side, we get 



1 dV,{p)dp{x,t) 
^p{x, t) J 2t dp dx 
(51) 

Considering Vc{p) — v'^ + Tf{p) and defining the "traffic 
pressure" as 

P{x,t)=p{x,t)e{x,t) + ^ (52) 
the corresponding macroscopic velocity equation becomes 



dt 
dp{p) d 



dx 



[p{x,t)V{x,t)\ -V{x,t) 



dp dx 

dp{p) dV{x,t) 



dp{p) dpjx, t) 
dp dx 



dp 



dx 



(57) 



By comparison with the macroscopic velocity equation 
((53)) we see that the model by Aw and Rascle does not 
have a relaxation term \Vc{p) — T^(a;,t)]/r, which would 
correspond to the limit t ^ oo. Moreover, we find 



idP{x,t) ^^dp{p) dv{x,t) 



p dx 



dp 



dx 



(58) 



Therefore, the traffic pressure according to the model of 
Aw and Rascle is a function of the velocity gradient rather 
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than the density gradient, in contrast to Payne's pressure 
term ([54| . Consequently, Aw's and Rascle's pressure term 
must result in a different way than Payne's one. In order to 
demonstrate this, let us now discuss a generalization of the 
optimal velocity model and its macroscopic counterpart. 



4.3.2 Non-Local Macroscopic Traffic IVIodels 

It is well-known that the optimal velocity model may pro- 
duce accidents, if the initial condition, the optimal velocity 
function Vo{d), and the parameter t are not carefully cho- 
sen. In order to have both, the emergence of traffic jams 
and the avoidance of accidents, we need to assume that the 
repulsive interaction force among vehicles does not only 
depend on the vehicle distance di{t) = Xi-i{t) — Xi{t), but 
also on the vehicle velocity Vi (t) (to reflect the dependence 
of the safe distance on the vehicle speed) or on the relative 
velocity 

Av.it) ^ v,{t) ^ v,-,it) ^ . (59) 

at 

The corresponding generalization of the acceleration equa- 
tion (fTTj) reads 



V^{t) 



f{d,{t),v,{t),Av,it)) . (60) 



This also changes the associated macroscopic traffic equa- 
tion. Namely, equation (|T7)) ') has to be replaced by 



dt 
1 



dx 



d 



p{x, t) dx 

2^ {p{x,ty 



[p{x,t)9{x,t)\ 



V{x,t) 



V{x,t),AV{x,t) 



Iff I 

2' \p{x + l/p,t)- 



V{x + l/p, t),AV{x + l/p,t) 



(61) 



For the sake of comparison with other macroscopic traf- 
fic models and linear stability analyses, let us perform a 
Taylor approximation of this. First, we may write 



1 



p{x + l/p,t) 
1 



,AV{x + l/p,t),V{x + l/p,t) 



p{x,t) 
df dd 
dd dp 
df 



AV{x,t),V{x,t) 



dv 
df 



dAv 



p{x + l/p,t) - p{x,t) 
V{x + l/p,t)~V{x,t) 
AV{x + l/p,t) ~ AV{x,t) 



Then, we may insert dd/dp = — l/p^, 
p{x + l/p, t) - p{x,t) 



dp_l 
dx p 



(62) 



(63) 



dV 1 

V{x + llp,t)-V{x,t)^—-. (64) 

ox p 

Furthermore, considering Z\ui(t) = —ddi/dt, p{x,t) — l/di{t), 
and the continuity equation dp/dt = dp/dt+ Vdp/dx = 
—pdV/dx, we get 



AVix^t)^-^^ ^ 



1 dp{x,t) 



and 



dt \p{x,t) J p{x,tY dt 

1 dV{x,t) 
p{x,t) dx 
V{x,t)-V{x + llp, t) 



AV{x + l/p, t)~ AV{x,t) 
dAV 1 1 d f dV 1 



(65) 



dx p p dx \ dx p ^ 
1 dpdV 1 d^V _ 1 d^V 
p^ dx dx p^ dx^ p^ dx^ ' 



(66) 



as a linearization drops products of gradient terms such as 
{dp/ dx){dV / dx) (which are assumed to be smaller than 
the linear terms). Altogether, with dd/dp — —1/p^ we can 
write 



1 



p{x + l/p, t) 
1 



.V{x + l/p, t),AV{x + l/p, t) 



I dfdV 1 df d'^V 



p dv dx p^ dAv dx^ 
With the definition 

Vo{p,V,AV)^v"+Tf (^^,V,AV 



(67) 



(68) 



we may finally write 

dV{x,t) , ,, dV{x,t) Id. . 



v,{p,v,Av)-v{x,t) 1 dfdp 

T 2p^ dd dx 

1 dfdV 1 df d^V 
2p dv dx 2p2 dAv dx"^ 



(69) 



Furthermore, let us assume that the variance can be ap- 
proximated as a function of the density and the average 
velocity: 

9ix,t)^9,{p{x,t),V{x,t)). (70) 
With the definitions 

dPi ...^ de,(p,V) , 1 df{\/p,V,AV) 

— =9.{p^V)+p^^ + —^ , (71) 



dP2^ de,{p,v) I df{i/p,v,Av) 

dV ^ dV 2 dv 
1 df{l/p,V,AV) 



2p2 



dAi 



(72) 
(73) 
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(where r/ should be greater than zero), we may also write 
the linearized macroscopic traffic equations as 



dV{x,t) 
dt 

_ldl\dp_ _ 

p dp dx 
Voip, AV,V) 



V{x,t) 



dV{x,t) 



dx 
1 dP2 dV 
p dV dx 
-V{x,t) 



' dx'^ 



(74) 



The term rjd^V/dx^ can be interpreted as viscosity term 
and has some smoothing effect. Further viscosity (and dif- 
fusion) terms may be derived by second-order Taylor ex- 
pansions. 

Note that the pressure term P2 looks similar to Eq. 
Therefore, it should be possible to derive a pres- 
sure term corresponding Aw's and Rascle's model from a 
suitable microscopic traffic model, but one would expect 
additional terms such as (|52p to occur as well. 



5 Summary, Discussion, and Conclusions 

In this paper, we have discussed several approaches to de- 
rive macroscopic traffic equations from microscopic car- 
following models. It has been pointed out that a Tay- 
lor approximation should be used only for linear stability 
analyses, as the gradients may otherwise be too large for 
the approximation to work. Further undesireable conse- 
quence of a gradient expansion are the possible occurence 
of negative velocities, diffusion instabilities, and inefficient 
numerical solution methods. 

The linear interpolation approach often works well in 
practise, but it is theoretically inconsistent as it violates 
the continuity equation which is required for the conserva- 
tion of the vehicle number. In constrast, the smooth par- 
ticle hydrodynamics approach was suited in all respects. 
It led to a non-local macroscopic traffic model, as did the 
gas-kinetic based traffic model. In order to have a realistic 
traffic dynamics (in particular accident avoidance if a ve- 
hicle with speed w*^ approaches a standing car), one needs 
to take into account that the repulsive vehicle interac- 
tions depend not only on the vehicle distance, but also on 
the relative velocity and the vehicle velocity. This leads 
to a specification of the traffic pressure which contains 
variance-dependent terms, additional terms proportional 
to dp/dx as in Payne's model, and further terms propor- 
tional to dV/dx as in Aw's and Rascle's model. While 
the variance-dependent term describes dispersion effects, 
Payne's, Aw's and Rascle's terms refiect effects of vehicle 
interactions. Finally note that, in case of multi-lane traffic, 
the additional inter-lane variance 



0{x,t) 



1 



_ Pi{x,t) 
p{x,t) 



Vi[x,t)-V{x,t)f 



(75) 



must be added to the inner-lane variance 9{x,t), where 
pi{x,t) is the density and Vi(x,t) the average velocity in 
lane I at location x and time t [2ll21l. 
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Christian Schmeiser at the Wolfgang Pauli Institute in Vienna 
from May 5-9, 2008, with partial support by the CNRS. 



A The Non-Local, Gas-Kinetic Based Traffic 
Model 

For comparison, let us shortly recall the form of the non- 
local gas- kinetic based traffic model (GKT model). This 
has been derived via a collision approximation, see Ref. 
[S] . This can be written in the form of equation ([55)1 with 
P{x,t) = p{x,t)8{x,t), but Vc{p) must be replaced by a 
non-local expression 



v^{p,v,9,p^,v^,e^ 



-t[i-p(pJ]x(p+)p+b(^) 



repulsive interaction term 



(76) 

Here, the index indicates evaluation at the advanced 
"interaction point" x + sq + TV, where Sq represents the 
minimum vehicle distance and TV the velocity-dependent 
safety distance. The related non-locality has some effects, 
which other macroscopic models generate by their pres- 
sure and viscosity terms. The dependence of the non-local 
repulsive interaction on the effective dimensionless veloc- 
ity difference 



A = 



V -V. 



6 - 2r, 



(77) 



takes into account effects of the velocity variances 6, 9^ , 
and velocity correlations r among successive cars 1211 . Fur- 
thermore, the "Boltzmann factor" 



B{A) = {9-2r 



AN{A) + {l + A^)E{A) 
(78) 

in the braking term is monotonically increasing with AV. 
It contains the normal distribution 



N{A) = 



V2^ 



and the Gaussian error function 

A 

E{A) ^ I dzN{z). 



(79) 



(80) 



To close the system of equations, the velocity correlation 
r is specified as a function of the density in accordance 
with empirical observations. Moreover, for a description 
of the presently known properties of traffic flows it seems 
sufficient to set 

9 = A{p)V^. (81) 

This guarantees that the velocity variance will vanish when- 
ever the average velocity goes to zero, but it will be posi- 
tive otherwise. It should be noted that the variance pref- 
actor A is higher in congested traffic than in free traffic. 
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The "effective cross section" is, finally, specified via 

[1 - p{p)]x{p) = ^ , (82) 

where T is the safe time headway and pj am the maximuni 
vehicle density. This formula makes also sense in the low- 
density limit p — > 0, where x 1 ^-nd p — s- 1. 

A linear stability analysis of the non-local traffic model 
can be done via a gradient expansion. It results in equa- 
tions of the kind ([74| and further viscosity and diffusion 
terms [22] . 
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